Static disorder and structural correlations in the low temperature phase of lithium 

imide 
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Based on ab-initio molecular dynamics simulations, we investigate the low temperature crystal 
structure of Li2NH which in spite of its great interest as H-storage material is still matter of debate. 
The dynamical simulations reveal a precise correlation in the fractional occupation of Li sites which 
leads average atomic positions in excellent agreement with diffraction data and solves inconsistencies 
of previous proposals. 



Lithium amide (LiNH 2 ) and imide (Li 2 NH) have been 
extensively studied in recent years as promising materi- 
als for hydrogen storage [l|-[5[. Hydrogen release occurs 
in the mixture LiNH2/LiH via a reversible solid state 
decomposition reaction into lithium imide and molecu- 
lar hydrogen, LiNH 2 + LiH -> Li 2 NH + H 2 . The typi- 
cal operating temperature for this system is around 280 
°C, which is probably too high for on-board applica- 
tions. Nevertheless, the amidc/imidc system is under 
deep scrutiny since it represents a prototypical, relatively 
simple system, which could shed light on the mecha- 
nisms of reversible H-rclease in the more complex, and 
technologically more promising, reactive hydrides made 
of mixtures of amide, borohidrides and/or alanates (e.g. 
LiNH 2 /LiBH 4 or LiNH 2 /NaAlH 4 )0, H- The search for 
better performing materials in this class would greatly 
benefit from a microscopic knowledge of the decompo- 
sition mechanism which requires in turn a detailed de- 
scription of the crystalline phases involved. For several 
complex hydrides the structure of the phases undergo- 
ing the dehydrogenation/rehydrogenation process is still 
not fully resolved. This is the case for instance of the 
most studied sodium alanate (NaAlB^) for which Raman 
spectroscopic data@ and ab-initio simulations^ very re- 
cently suggested the existence of a new high temperature 
phase which is expected to mediate the decomposition re- 
action in place of the low temperature a-phasc considered 
so far. The structure of Li 2 NH as well is still a matter of 
debate. Structural refinement from neutron and x-rays 
diffraction data reveals a structure of Li 2 NH with frac- 
tional occupation. In spite of a substantial amount of ex- 
perimental and theoretical investigation, the problem of 
the actual local structure which yields this long-range dis- 
order is still unsettled. Based on ab-initio simulations, we 
have identified a model for the local structure of the low- 
temperature phase of lithium imide (Li 2 NH) which solves 
inconsistencies of previous proposals and fully agree with 
experimental data available. 

Differential thermal analysis and NMR 
measurements [1, Q in the late 60's revealed a re- 
versible phase transition at 356 K between an unknown 
low temperature (LT) structure and a high tempera- 
ture antifluorite phase of Li 2 NH. Structure refinement 



from x-ray and neutron diffraction measurements on 
deuterated imide (Li 2 ND) have been published only 
very recently flfjj. The high temperature phase yields 
a diffraction pattern consistent with an anti-fluorite 
structure, in which hydrogen atoms occupy the 192Z 
positions of the Fm3m space group. At low temperature 
(100-300 K) the diffraction data were best fitted by a 
cubic crystal with Fd3m space group. The LT crystal 
was described as a superstructure of the antifluorite 
phase in which one out of 8 Li atoms is displaced to 
an interstitial site giving rise to Li vacancies arranged 
in an ordered manner and tetrahedrally coordinated to 
four NH groups (Fig. [Ha)). The structure is stabilized 
by electrostatic interaction of H 5+ pointing towards the 
Li + vacancy which is formally a negatively charged site. 
Similar tctrahedral arrangements of four NH groups are 
present also in other imides, such as Li 2 Mg(NH) 2 [llj. In 




Figure 1. (a) Tetrahedral arrangement of NH groups around 
a Li + vacancy in the LT phase of Li2NH. (b) Li ions are dis- 
tributed on the three sites independent by symmetry labeled 
in panel. The interstitial Li3 sites have fractional occupation. 
Sites which are left empty in the Ima2 model of Ref. [I^l are 
shown as transparent spheres, (c) Occupied Li4 octahedral 
sites obtained upon relaxation of the Ima2 model, (d) Our 
proposed correlation in the occupation of Li3 and Li2 sites 
leads to a tetrahedral cluster of four interstitials ions at Li3 
sites around a vacancy at Li2 site. 
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the LT phase, H and N atoms occupy 32e sites. Li atoms 
are distributed over three different sites (Fig. [ljc)): 48/ 
(Lil), 8a (Li2) and 32e (Li3). The Li3 site corresponds 
to a displacement along (111) directions from the 
octahedral position Li4 (sites 16d of the Fdim space 
group, Fig. 1(c)). According to the Rietveld refinement, 
the latter has a fractional occupation of about 1/3 at 
low temperatures, which only at room temperature gets 
closer to the value required by stoichiometry (1/4). 

Diffraction data, however, do not provide informa- 
tion on correlations among the occupation of sites 32e 
which are mandatory to get full insight on the structural 
and dynamical properties of the system. To this aim, 
an attempt was made by Herbst and Hector to model 
the Fd2>m crystal by fully occupying selected 32e sites 
III 12, 111 ■ In the structure proposed, a single Li + ion 
is present around each Li2 site, resulting in a Ima2 sym- 
metry (Fig. [11(b)). The distance between Li2 and Li3 
sites is, however, very small (1.53 A) and upon geometry 
optimization performed by ab-initio calculations fl2T| . Li 
interstitials are displaced from Li3 sites into octahedral 
Li4 positions (sites 16c? of the FdZm space group) with a 
longer Li4-Li2 distance of 2.3 A (Fig. WLc))- The locally 
relaxed structure obtained by Hector and Herbst [12j 
has Imma symmetry and will be referred to hereafter 
as the HH phase. The positions of Li interstitials at 
Li4 sites are, however, incompatible with Rietveld re- 
finement. Moreover, the theoretical formation energy 
of this phase is too high when compared to experimen- 
tal data. Thus, other structures have been proposed on 
the basis of theoretical calculations, with Li atoms ar- 
ranged in an ordered manner and a slightly lower for- 
mation energy [li, 15 1. However, their equilibrium lattice 
parameters are inconsistent with the space group symme- 
try inferred experimentally. On the other hand, the local 
instability of the ideal antifluorite structure of Li ions 
at low temperature was confirmed by ab-initio molec- 
ular dynamics simulations at 300 K which showed the 
spontaneous formation of Li Frenkel pairs [l(| and strong 
distorsions of the Li sublattice 17\. 

Based on ab-initio simulations, we propose a structure 
for the LT phase which solves the problems mentioned 
above by introducing vacancies on partially occupied Li2 
sites with a precise correlation with the occupation of Li3 
interstitial sites. The solution of the puzzle came from 
the analysis of ab-initio molecular dynamics trajectories. 

We started our analysis by performing ab-initio molec- 
ular dynamics simulations on the HH structure. A s/2 x 
1 x v2, 128-atoms supercell corresponding to two Imma 
unit cells was used. We performed Born-Oppenheimer 
molecular dynamics simulations within Density Func- 
tional Theory (DFT) with gradient corrected exchange 
and correlation functional j 1 8| as implemented in the 
CPMD package. UltrasoftflJ and Goedecker- 
tyP e (m pseudopotentials were used, respectively, for N 
and H atoms and for Li with three valence electrons. 



Kohn-Sham orbitals were expanded in plane waves up 
to a kinetic-energy cutoff of 50 Ry. Brillouin Zone (BZ) 
integration was restricted to the L point only. A time 
step of 0.6 fs was used and a constant temperature of 
300K was enforced by an optimal-sampling generalized 
Langevin equation thermostat j22j . Equilibrium geome- 
tries of relevant structures emerged from the dynamical 
simulations were optimized with special k-points meshes 
and the Quantum-Espresso suite of programs (2^1 . Ac- 
tivation energies for diffusion processes discussed below 
were obtained by Nudged Elastic Band optimizations [24| . 

In the dynamical simulations, 15 ps long, we observed 
that constitutional vacancies coordinated by NH groups 
are very stable, the latter performing only small libra- 
tions around their equilibrium position. Instead, Li in- 
terstitials occupying Li4 positions are very mobile and 
we observed several jumps between adjacent sites, tak- 
ing place via an exchange mechanism with one of the 
Lil atoms (cf. Fig. SI in supplementary materials be- 
low). The activation energy for diffusion of interstitials 
turns out to be as low as 0.13 eV. Moreover, we verified 
that many different arrangements of interstitial atoms 
in Li4 sites have total energies within few meV. At finite 
temperature, we would thus expect a disordered arrange- 
ment of interstitial atoms in the Li4 sites, leading to a 
fractional occupation of 1/2 and an overall Fd3m space 
group. 




Figure 2. a) The 192 atoms supercell is shown (3x1x1 HH 
unit cells). The labeling indicates all possible sites where the 
four Li2 vacancy - Li3 interstitials clusters can be distributed 
inside the supercell. We found that the lowest-energy arrange- 
ments are ACFL (-30 meV/f.u., see also Fig. S2 below) and 
AEGK (-29). Energies are given relative to the relaxed HH 
structure. Configurations where the tetrahedra are unevenly 
distributed such as AEIC (-17) or ACEG (4) are higher in 
energy, and those where nearest-neighbor sites are occupied 
simultaneously are even less stable. 

However, the positions of Li4 still differ from those of 
Li3 obtained from Rietveld refinement. Moreover, the 
low diffusion barriers found would assign a supcrionic 
character to Li2NH even at room temperature for which 
no experimental evidence has been provided so far. 

The solution to these discrepancies came from a closer 
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inspection of the molecular dynamics trajectory, which 
revealed that as soon as three interstitial atoms happen 
to be all nearest neighbors to an occupied Li2 site, the 
latter atom is displaced from its position. A tetrahc- 
dral cluster is formed with a vacant Li2 site at its cen- 
ter (Fig. QJd)). These clusters remain stable thereafter 
and prevent diffusion of Li interstitials. The interstitial 
Li now occupies the Li3 32e sites and a vacancy is in- 
troduced at site Li2 in such a way that the very short 
Li3-Li2 distance of the unrelaxed HH model is removed. 

Having observed the propensity of Li atoms to form 
these rather stable tetrahedral structure, we propose a 
new stable structure in excellent agreement with diffrac- 
tion data, as follows. At T=0 we expect that all the 
Li interstitials are involved in the formation of tetrahe- 
dral clusters. The smallest supercell compatible with 
this requirement and consistent with the stoichiometry 
is formed by 3x1x1 HH unit cells (192 atoms). In such 
supercell there will be four vacancy-interstitials clusters, 
which can be distributed over twelve Li2 sites (see Fig- 
ure [2]) . We considered all the possible arrangements of 
the vacancy-interstitials tetrahedra independent by sym- 
metry excluding those where two first-neighbor Li3 sites 
are occupied. 

Provided that the clusters are evenly distributed in the 
cell to balance the negatively-charged constitutional va- 
cancies, the spread in the total energy of the different con- 
figurations is well within thermal energy at room temper- 
ature (see Fig. We repeated our analysis on a larger 
3x1x3 supercell finding similar results. Such degener- 
acy allows for long-range disorder in the arrangement of 
tetrahedra, so that a diffraction pattern consistent with 
a higher-symmetry Fd3m, with a fractional occupation 
of sites Li3 (1/3) and Li2 (2/3) is to be expected. These 
fractional occupations are compatible with the stoichiom- 
etry and solve the inconsistency in the model of Ref. [l(J ■ 

By averaging over possible configurations of our 192- 
atoms supercell we obtained the symmetry-adapted av- 
erage positions reported in Table U and the relative mean 
square displacements which would correspond to a static 
contribution to the Debye- Waller factor. The presence of 
clusters generates long-range distortions in the Lil lattice 
and in the orientation of H atoms around a constitutional 
vacancy, resulting in a large mean square displacement 
from the average position for these species. The thermal 
Debye- Waller factors as a function of temperature have 
been computed from harmonic phonons (at the super- 
cell T-point) and are given in Fig. S3 in supplementary 
materials below and at T = 100 K in Table Q] 

The agreement with the experimental positions from 
Rietveld refinement is excellent. The diffraction pattern 
of our proposed structure compares well with that de- 
duced from the experimental positions and fractional oc- 
cupation (Fig. S4 in supplementary materials below). 

The formation enthalpy of the ACFL configuration 
(see Fig. [5]) is lower (-30 meV f.u. -1 ) than that of the 



Site 



(exp. Ref [10]) Ui SO static/thermal 



N 32e f=l 
D 32e f=l 
Lil 48/ f=l 
Li2 8a f=2/3 
Li3 32e f=l/3 



0.2408 (0.2418) 
0.2982 (0.2976) 
0.3755 (0.3734) 

0.0367 (0.0376) 



0.34/0.50 
1.58/2.00 
2.18/0.93 
0.25/0.87 
0.22/0.91 



Table I. Structural parameters of the best-fit Fd3m structure 
corresponding to the ACFL arrangement of tetrahedral clus- 
ters in a 3 x 1 x 1 supercell (see Fig. [2]). The Debye- Waller 
factor (1CT 2 A 2 ) is the sum of a contribution due to static dis- 
order, and a thermal contribution calculated at T = 100 K 
for Li2ND. The experimental data from Rietveld refinement 
of neutron diffraction patterns at 100 K from Ref. [l^l are 
given in parentheses. 



HH structure. Considering also zero-point energy and 
the harmonic contribution to the total free energy, our 
structure is still more stable than HH's (-18 meV f.u. -1 
at 298 K) and the structure proposed in Ref. [3] (-3 meV 
f.u. -1 ). Our structure is slightly less stable (by 22 meV 
f.u. -1 ) than that suggested in ref.|15( which however does 
not correspond to the experimental space group. How- 
ever, uncertainties in the relative energies of the different 
phases are expected due to current approximations to 
the exchange and correlation functional as occurs for in- 
stance in other materials we have recently studied pB|. 
Moreover, entropic contributions due to disorder might 
compete with small enthalpy differences of ordered mod- 
els and need to be taken into account before assigning the 
most stable phase at the DFT level. Because of all the 
above limitations, our aim was not to find the lowest pos- 
sible ordered structure on the potential energy surface as 
pursued in Refs. [H, 15 but to refine a disordered struc- 
ture which was a local minimum in DFT in such a way 
that is compatible with the experiments. 

We verified the stability of the tetrahedral clusters by 
performing 20 ps of molecular dynamics at 300 K us- 
ing the most stable arrangement ACFL. The vacancy- 
interstitials clusters do not break nor constitutional va- 
cancies diffuse on the simulation time scale. The com- 
puted energy barrier to break a tetrahedral cluster by 
diffusion of a Li3 ion is indeed 0.5 eV, as obtained by 
NEB optimization, resulting in a low mobility of Li atoms 
in the LT phase (Li3 diffuses by an exchange mechanism 
similar to that of Li4, see Fig. SI in supplementary ma- 
terials below). The energy cost to break a cluster in the 
ACFL structure, via the reaction 4Li3 — > Li2 + 3Li4 is 
AE=0.20 eV. This energy depends on the final configu- 
ration of the interstitials, and on the number of clusters 
which arc simultaneously broken. For instance, the dif- 
ference in energy with the HH structure - with all the 
clusters broken - amounts to AE=0.36 eV per cluster. 

We can compute the change in configurational en- 
tropy due to breaking of tetrahedra and formation of 
these interstitial Li4 atoms, by assuming a perfectly ran- 
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dom distribution of tetrahedra and occupied Li4 inter- 
stitial sites (see supplementary materials below) which 
finally yields a concentration of broken tetrahedra of 
x ~ \/4/3 exp[— AEo/3fcsT]. This corresponds to a con- 
centration at 300 K between ~6% and ~0.5%, depend- 
ing on the estimate taken for AEo. In spite of this large 
concentration of Li4 the diffusivity is still low because 
most of the percolating paths for Li4 interstitials are sup- 
pressed by the presence of Li3 atoms fixed in tetrahedra. 

We speculate that the possible coexistence of partially 
occupied Li3 and Li4 sites might contribute to the change 
with temperature of the occupation of Li3 sites emerged 
from the Rietveld analysis • We encourage new refine- 
ment of the diffraction data starting from our structural 
parameters in Table HI eventually allowing for a (small) 
partial occupation of sites Li4. 

In summary, by combining ab-initio molecular dy- 
namics and structural optimization, we have provided 
a full description of the structure of the LT phase of 
Li2NH. Atomic positions coincide with those inferred 
from diffraction data but partial occupations of Li sites 
have a strong nearest neighbor correlation which solve 
the inconsistencies raised by previous proposals. This 
work represents an exemplary demonstration of how dy- 
namical simulations can provide crucial insights to fully 
resolve the structure of systems with partial disorder, 
complementing experimental diffraction data. 

We gratefully thank W. I. F. David for discussion and 
information. 
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Supplementary materials 



COMPUTATIONAL DETAILS 



We started our analysis by performing ab-initio molecular dynamics simulations on the HH structure. A y/2x 1 x s/2, 
128-atoms supercell corresponding to two Imma unit cells was used. We performed Born-Oppenheimer molec- 
ular dynamics simulations within Density Functional Theory with gradient corrected exchange and correlation 
functional [HI as implemented in the CPMD [n| package. Ultrasoft[2(| and Goedecker-typejilj pseudopotentials 
were used,respectively, for N and H atoms and for Li with three valence electrons. Kohn-Sham orbitals were ex- 
panded in plane waves up to a kinetic-energy cutoff of 50 Ry. Brillouin Zone (BZ) integration was restricted to the T 
point only. A time step of 0.6 fs was used and a constant temperature of 300K was enforced by an optimal-sampling 
generalized Langevin equation thermostat 22j. Equilibrium geometries of relevant structures emerged from the dy- 



namical simulations were optimized with special k-points meshes and the Quantum-Espresso suite of programs 
Activation energies for diffusion processes discussed below were obtained by Nudged Elastic Band optimizations 



23] 
24] 



MECHANISM OF DIFFUSION OF LI INTERSTITIALS 




Figure SI. Diffusion exchange mechanism for Li interstitial. In the picture the initial, transition and final states are shown. 
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DEBYE- WALLER FACTOR 



We calculated the Debye- Waller factor for each species as a function of the temperature from harmonic phonons. 
The mean square displacement as a function of the temperature in harmonic approximation is given by 

2 1 ^ H |e(m,z)| 2 [ / hu; m \ ll 

m,i 

were a runs over the atomic species, M a is the mass of ath species, i runs over A^ a atoms of species a, while u) m and 
e(m, i) are the frequency and the eigenvector of the mth harmonic phonon. 




Figure S2. Debye-Waller factor, U iso = l/3{u 2 a ), for D, N, Lil, Li2, Li3 as a function of the temperature. 
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LI 2 NH X-RAY DIFFRACTION PATTERN 



We show in Fig. n a comparison of the calculated x-ray diffraction pattern using both experimental [10( and our 
theoretical positions. 



0.05 i 

: 

-0.05 ; 

1 - 



r— 



th. 

exp. 




I. 



L 



50 



150 



Figure S3. Bragg reflections intensities computed the average positions and occupations reported in Table 1, compared with 
those obtained using experimental values. No instrument-related corrections have been applied and the same Debye- Waller 
factor (Uiso = 0.05 A 2 ) has been used for all atomic species. 



CONFIGURATIONAL ENTROPY OF TETRAHEDRAL CLUSTERS 



To compute the change in computational entropy due to the breaking of tctrahcdra we considered that the number 



of tetrahedra Nt is related to the number of Li4 and Li2 sites by Nm 
energy due to the breaking of a fraction x of tetrahedra is given by 
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which yields an equilibrium concentration of 



exp[— AE /3fesT]. 



